Core reconstruction of the 90° partial dislocation in non-polar semiconductors 
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We investigate the energetics of the single-period and double-period core reconstructions of the 90° 
partial dislocation in the homopolar semiconductors C, Si, and Ge. The double-period geometry 
is found to be lower in energy in all three materials, and the energy difference between the two 
geometries is shown to follow the same trends as the energy gap and the stiffness. Both structures 
are fully reconstructed, consisting entirely of fourfold coordinated atoms. They differ primarily 
in the detail of the local strains introduced by the the two reconstructions in the core region. The 
double-period structure is shown to introduce smaller average bond-length deviations, at the expense 
of slightly larger average bond-angle bending distortions, with respect to the single-period core. 
The balance between these two strain components leads to the lower energy of the double-period 
reconstruction. 
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A fundamental understanding of plasticity in solids 
clearly requires a knowledge of the atomistic structure at 
the cores of dislocations. In particular, the microscopic 
mechanisms of dislocation motion are intimately related 
to the defects (e.g., kinks) that can occur within the dis- 
location, which in turn are connected with the underlying 
lattice symmetries and the nature of the reconstruction 
in the core.u Recently, advances in computer power and 
computational methodology have led to an active area of 
research focused on the theoretical study of the atomistic, 
structure of these dislocation cores and their defects.BrEil 
However, in at least one important case, even the struc- 
ture of the dislocation core itself remains fundamentally 
in doubt. 

In tetrahedrally bonded semiconductors, the two most 
frequently occurring dislocations are the 30° and the 90° 
partial dislocations, lying on {111} planes along [110] 
directions .EjCj These materials are of obvious techno- 
logical importance, and the detailed understanding of 
the atomic structure at the dislocation cores is of great 
interest, since dislocations influence both the electronic 
and the mechanical properties of semiconductor devices. 
Indeed, a great deal of theoretical effort has been de- 
voted to fttudy the properties of dislocations in these 
materials.ErEiJ In the particular case of silicon, the the- 
oretical study of the dislocation cores at the atomistic 
scale has revealed a rich structure of point excitations 
(kinks and recon&tpniitinn defects) in the core of the 30° 
and 90° partials.p~BE3T3 

Most of these works have concentrated on understand- 
ing the structure of the 90° partial dislocation. Un- 
til recently, a consensus had emerged at the theoreti- 
cal level, about the nature of the reconstruction at the 
core of the defect. In its unreconstructed configuration, 
the core of the 90° partial displays a zigzag chain of 
threefold-coordinated atoms running along the disloca- 



tion direction, with broken bonds lying nearly parallel 
to the slip plane. Mirror symmetry planes are present is 
this configuration, as can be seen in Fig. 0(a). A vari- 
ant of this structure is one in which the dashed lines in 
Fig. 0(a) are considered to be covalent bonds, resulting in 
a "quasifivcfold" reconstruction that also retains the mir- 
ror symmetry On the other hand, a reconstruction that 
breaks the mirror symmetry of the unreconstructed core, 
while preserving the lattice periodicity along the line, is 
shown in Fig. 0(b). In this case all dangling bonds have 
been eliminated, and all the atoms are fourfold coordi- 
nated. Such a reconstruction was predicted to be sub- 
stantially lower in energy than both the unreconstructed 
and the quasifivefold-reconstructed cores j3i3'Eirt2lll3 and 
thus to be the one expected to occur in nature. Wc 
will refer to this symmetry-breaking reconstruction as the 
single-period (SP) reconstruction. 

However, our recent theoretical work on the core re- 
construction of this dislocation in Si produced a surprise. 
In Ref. ||, we proposed an alternative core structure for 
the 90° partial in which, in addition to symmetry break- 
ing of the SP core, the periodicity along the dislocation 
line is doubled. This double-period (DP) reconstruction, 
which is shown in Fig. |l|(c), can be derived from the 
SP one by introducing alternating kinks at every lattice 
site along the core. This geometry is consistent with all 
available experimental information about the 90° par- 
tial. Like the SP core, the DP structure is fully re- 
constructed, and thus neither one gives rise to deep-gap 
states which would show an EPR signal. EPR exper- 
iments in Si indicate a very smalL density of dangling 
bonds in the core of dislocation]^ Moreover, both 
cores consist entirely of fivefold, sixfold, and sevenfold 
rings, both being consistent with images produced by 
transmission electron microscopy, at the current level of 
resolution of this technique.^ The results we obtained in 
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FIG. 1. (a) Symmetric reconstruction of the 90° partial 
dislocation in homopolar semiconductors. Shaded area in- 
dicates stacking fault, (b) The single-period (SP) symme- 
try-breaking reconstruction. (c) The double-period (DP) 
symmetry- breaking reconstruction. 

Rcf. H show the DP structure to be lower in energy than 
the SP one, by means of Keating-potential, total-energy 
tight-binding (TETB), and ab initio local-density (LDA) 
calculations. 

In the present work, we investigate the issue of DP 
versus SP reconstruction in the homopolar diamond- 
structure semiconductors C, Si, and Ge. As in our pre- 
vious work on Si, we present LDA, TETB, and Keating- 
potential results for the energies of the two core recon- 
structions. The trend in the energy difference between 
the DP and SP cores is seen to correlate with the stiff- 
ness of each material, with the DP core being more (less) 
strongly favored in C (Ge), as revealed by our ab-initio 



TABLE I. Calculated energy difference in meV/A, between 
the SP- and DP-core reconstructions of the 90° partial in C, 
Si, and Ge. Cell size refers to the double-period cell. £bp 
is the energy of the double-period reconstruction. For the 
single-period case, _Bsp and AEsp are respectively the average 
and difference of the energies for the two different relative 
arrangements of mirror symmetry-breaking. 
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34 
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a Evaluated at LDA-relaxed structure. 



and TETB calculations. This is in contrast with the 
Keating-potential calculations which, despite predicting 
the correct ground-state and trend between Si and Ge, 
favor the SP structure over the DP core in the case of C. 

As in Ref. [|, we employed supercells of 96 and 192 
atoms for the SP and DP cores respectively, containing 
a dislocation dipolc in the quadrupole arrangement sug- 
gested in Ref. [| These were fully relaxed (with average 
forces no larger than 0.01 eV/A in each case) using an 
LDA approach, with core states represented by Troullier- 
Martins pseudopotentials, as implemented in the fhi96md 
package. lj All the energy differences were converged to 
within ~5% with respect to plane-wave cutoff. For Si 
and C, these cells were also relaxed with a-TETB model, 
using an 0(N) density-matrix techniqueEZl to solve for 
the electronic structure. (The corresponding TETB re- 
sults for Ge would be expected to resemble tbose for Si.) 
The tight-binding Hamiltonian of Kwon et al£B was used 
for Si, while the model proposed by Xu et al. was ap- 
plied to Cc3 The convergence of the TETB results with 
respect to cell size was checked by relaxing larger cells, 
with 576 (288) atoms for the DP (SP) structure. Keating- 
potential results were used in order to investigate qualita- 
tively the strains associated with the two structures. For 
these, we-use the original set of parameters introduced by 
KeatingcJ for the three materials. Further details about 
the supercells and the technical aspects of our TB calcu- 
lations can be found in Refs. |^-|4. 

We discuss first the LDA results for the 192-atom su- 
percell as shown in Table I. The calculated energy of 
the SP core depends on whether the breaking of mirror 
symmetry occurs with the same sense, or with opposite 
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sense, for the two dislocations in the supercell. (There 
is a corresponding energy splitting in the DP core, but 
this is of a much smaller magnitude and was not taken 
into consideration.) As in Ref. |^, the average energy of 
these two possibilities is denoted by Egp, while the dif- 
ference is denoted by AEgp; the latter is expected to 
vanish in the limit that the supercell gets large. Focus- 
ing for the moment on the Egp values, it can be seen that 
the DP structure is preferred by a substantial margin for 
all three, and that this preference follows the same trends 
as the stiffness and the size of the band gap of the ma- 
terial. That is, C shows the strongest tendency toward 
stabilization of the DP core, followed by Si and then Ge. 

Since we cannot easily afford to repeat the LDA calcu- 
lation for a larger supercell, we have carried out parallel 
calculations on 192-atom and 588-atom supercells of C 
and Si using the TETB method. Table I shows first of 
all that the TETB results are in qualitative agreement 
with the LDA ones for both materials, although in the 
case of C the corresponding values are underestimated 
by a factor of about two. Secondly, the TETB results 
give us a good estimate of the importance of the finite 
supercell-size effects. That is, by inspecting the TETB 
results for the two supercells, it appears that one can say 
that E$p gives the energy of the true isolated SP struc- 
ture to within an error bar of approximately ±A£?sp/2. 
Applying this heuristic to the LDA results, we see that 
the supercell-size error is almost certainly insufficient to 
reverse the sign of the predicted £"dp — Egp- Returning 
now to the discrepancy between the TETB and LDA re- 
sults for C, this could be attributed at first sight to the 
fact that the model of Xu et al. that we use for C re- 
sults in elastic constants which are too soft [{Cxi — C12) 
and C44 are smaller by 35% and 17% with respect to ex- 
perimental values, respectively]. On the other hand, the 
TETB results for Si that we obtained using the model of 
Ref. [31], which also underestimates the elastic constants, 
are practically identical to those shown in Table I for the 
Kwon model. 

In order to gain some further insight on the difference 
between the SP and DP structures, we look at the maxi- 
mum deviations of bond lengths and bond angles for the 
LDA geometries, as shown in Table II. The main trend 
as a function of the material considered is the smaller 
variation of bond angles, and greater variation of bond 
lengths, in C relative to Si and Ge. This is consistent with 
the fact that the ratio of bond-angle bending to bond- 
stretching forces is bigger in C, compared to Si and Ge. 
However, we are mainly interested in the trends in going 
from the SP to the DP structure. For all three materials, 
we find that the pattern of deviations looks surprisingly 
similar for the two structures. In particular, the maxi- 
mum and minimum bond-length and bond-angle varia- 
tions are not systematically smaller in the DP structure, 
in spite of its lower energy. However, a trend becomes 
more visible when inspecting the root-mean-square devi- 
ations: we see that the rms bond-length variations are 
systematically slightly larger in the SP core. While it 



TABLE II. Minimum, maximum, and root-mean-square 
variations of bond lengths and bond angles for the 
LDA-relaxed SP and DP structures, relative to the corre- 
sponding bulk diamond values. 
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+3.1% 


+3.5% 
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2.8% 


2.5% 


4.0% 


4.1% 



also appears that the bond-angle variations go the other 
way, being slightly smaller for the SP core, this seems to 
be a weaker effect. These results suggest that advantage 
of the DP structure is that it allows a better packing of 
the atoms, as measured by bond-length variations, at the 
expense of a slightly greater bond-bending strain. 

In view of the fact that the variations in bond lengths 
and bond angles between the SP and DP structures are 
so subtle, it is not at all clear whether a Keating model 
could be expected to reproduce the correct trends in E^p 
vs. Ep>p for the three materials. In order to test this, we 
also present in Table I the energies for the two core struc- 
tures as computed using the Keating potential. This was 
done both by consistently relaxing and evaluating the en- 
ergy using the Keating model, and by simply evaluating 
the Keating energy of the LDA-relaxed supercells. Note 
that in the case of C, the DP-SP energy difference for the 
Keating-relaxed geometry is substantially smaller than 
the corresponding LDA result. Also, this set of numbers 
predicts an incorrect trend for this quantity among the 
three materials, as compared to the LDA results. For 
the larger supercell this quantity is even in qualitative 
disagreement with the TETB result, in the case of C. 
(We evaluated the Keating energy for the TETB geome- 
tries in these larger cells, and found the correct qualita- 
tive behavior, with Ep,p — Egp = —85 meV/A for C.) 
On the other hand, for the LDA geometries, the Keating 
results follow the same trends as those obtained with the 
other methods. From these results, we conclude that the 
Keating potential cannot be trusted to capture the en- 
ergy difference between the SP and DP structures on a 
quantitative (or perhaps even qualitative) level. 

To summarize, we investigated the energetics of the 
SP and DP core reconstructions of the 90° partial dislo- 
cation in the homopolar semiconductors, C, Si, and Ge. 
We find the DP core to be favored in all three materials, 
and observe that the energy difference between the two 
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geometries follows the same trends as the energy gap and 
the stiffness in these materials. The nature of this differ- 
ence is primarily associated with the local strains intro- 
duced by the the two reconstructions. The DP structure 
is shown to introduce smaller average bond-length devi- 
ations, at the expansion of slightly larger average bond- 
angle bending distortions, with respect to the SP one, 
with the delicate balance between these two strain com- 
ponents favoring the DP reconstruction. 
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